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The existence of extra degrees of freedom beyond the electroweak scale may allow the formation of 
black holes in nearly horizontal neutrino air showers. In this paper we examine the average properties 
of the light descendants of these black holes. Our analysis indicates that black hole decay gives rise to 
deeply penetrating showers with an electromagnetic component which differs substantially from that 
in conventional neutrino interactions, allowing a good characterization of the phenomenon against 
background. Naturally occurring black holes in horizontal neutrino showers could be detected and 
studied with the Auger air shower array. Since the expected black hole production rate at Auger is 
> 1 event/year, this cosmic ray observatory could be potentially powerful in probing models with 
extra dimensions and TeV-scale gravity. 



One of the most outstanding phenomena of TeV-scale 
quantum gravity Q is the possible production of semi- 
classical black holes (BHs) in particle collisions 
These BHs are expected to decay promptly, giving rise 
to large multiplicity events with large total transverse 
energy and a characteristic ratio of hadronic to leptonic 
activity of roughly 5:1. Production rates for the Large 
Hadron Collider (LHC) are found to be sizeable for a fun- 
damental Planck scale M* = 1 TeV ||||] . Additionally, 
BHs occurring very deep in the atmosphere (revealed as 
intermediate states of ultra high energy neutrino inter- 
actions) may trigger quasi-horizontal showers that could 
be detected with the Auger Observatory ||. The goal 
of this paper is to point out some salient experimental 
signatures of these air showers. 

We start the discussion by reviewing the relevant BH 
properties. BHs are believed to be described by semi- 
classical general relativity when their mass Mbh S> M». 
As Mbh approaches M», string excitations can become 
important and the BH properties rather complex. The 
ensuing discussion will be framed in the context of fiat 
extra dimensions, and we will comment on the warped 
scenario after presenting our results. In what follows, 
we rely on simple semiclassical arguments assuming that 
stringy effects are under control if Mbh/M* > 5. The 
Schwarzschild radius Rs of a (4+n) dimensional BH is 
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Hence, if one envisions a head-on collision involving par- 
tons i and j with cm. energy y§ = Mbh and impact 
parameter less than Rs, semiclassical reasoning suggests 
that a BH is formed. The total cross section of the pro- 
cess can be estimated from geometrical arguments 
and is of order 



(s) a nRl 



1 

M 2 



Mi 



BH 



8T(2i±3) 



M* 



(2) 



Before proceeding, we take note of a serious challenge to a 
geometric cross section raised by Voloshin |Q . The criti- 
cism centers on the exponential suppression of transitions 
involving a (few-particle) quantum state to a (many- 
particle) semi-classical state. In response ||, the geo- 
metric result was reaffirmed by arguing that it connects 
smoothly to the string scattering cross section in an en- 
ergy regime characterizing the transition to black hole 
physics. Whichever point of view one may find more 
convincing, it seems most conservative at this point to 
depend on experiment (if possible) to resolve the issue. 

With this in mind, the neutrino nucleon cross section 
reads Q 

® vN — >bh = Y2 I dx ai(xs) Mx,Q 2 ) , (3) 

where s = 2m^E Vl fi{x,Q 2 ) are parton distribution 
functions (PDFs), Mj^ 1 is the minimum BH mass, and 
the sum is carried out over all partons in the nucleon. 
Following m, the cross section is calculated using the 
CTEQ5MI PDFs §] with the momentum transfer Q 
taken to be equal to Mbh = \/xs- 

The energy released in neutrinos by supernova ex- 
plosions imposes several constraints on the fundamen- 
tal Planck scale |o|. Namely, M* > 500 - 1600 TeV, 
M* > 7-60 TeV and M* > 1 TeV, for n = 2, 3, 4, respec- 
tively. Therefore, a straightforward calculation shows 
that a v at^bh > &f,N ■> if > 4. Here, 
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is the total charged current Standard Model vN cross 
section (10 16 eV < E„ < 10 21 eV) @. For M BH « 5 
TeV, M* = 1 TeV, and neutrino primary energies around 
10 20 eV, c„jV-»BH ^ 10 31 cm 2 . Note that although the 
atmosphere presents a target of thickness of about 1000 
g/cm 2 to particles arriving vertically, the thickness in- 
creases up to ~ 36000 g/cm 2 to those arriving tangen- 
tially to the earth surface (i.e., with horizontal incidence 
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to the ground). Consequently, the probability of BH pro- 
duction is not negligible. Specifically, more than one BH 
event per year (2 events/yr for n = 4) could be detected 
by the ground array of the Auger Observatory 

The BH lifetime, governed by the Hawking evaporation 
process, is M 
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Thus, for Mbh > M*, t bh M B h > 1, the BH is a well de- 
fined resonance and may be thought as an intermediate 
state in the s— channel. Therefore, if one assumes that 
the BH evaporates instantaneously at its original temper- 
ature into its decay products, the average multiplicity is 
roughly M BH /(2T H ) [|, or equivalently, 



(N) 



1 



Mi 



BH 



■i + 2 



(6) 



where T B is the Hawking temperature. 

Most of the large multiplicity of observable quanta 
emitted in the BH decay is expected to come through 
hadronic jets produced by the quarks. The precise na- 
ture of the fragmentation process is unknown. We shall 
use here the quark — > hadron fragmentation spectrum 
originally suggested by Hill jjjl 
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that is consistent with the so-called "leading-log QCD" 
behavior and seems to reproduce quite well the multi- 
plicity growth as seen in colliders experiments. Here, x = 
E/Ej e t, E is the energy of any hadron in the jet, and E- }C t 
is the total energy in the jet. With the infrared cutoff set 
to x = 10 -3 , the average multiplicity per jet is approxi- 
mately 54 . The main features of the jet fragmentation 
process derived from dN^/dx s» (15/16) x~ 3 / 2 (1 — a;) 2 
(which provides a reasonable parametrization of Eq.(M) 
for 10~ 3 < x < 1) are listed in Table I. Now, assum- 
ing that the BH decays into all Standard Model particles 
(with equipartition among the particle species) and that 
each quark produces one hadronic jet, we obtain the "vis- 
ible" BH decay spectrum. 

We turn now to the analysis of the atmospheric cascade 
development triggered by the BH secondaries. In order 
to propagate the particles in the atmosphere we use the 
algorithms of AIRESQ (version 2.1.1) |l3]] . The showering 
of each charged hadron in the spectrum is simulated by a 
proton cascade of energy E, whereas the shower induced 
by a 7T° decay Q is replaced by a superposition of 2 pho- 
ton showers of energy E/2. In ~ 1/3 of the events, the 
lcptonic channel generates electromagnetic showers sim- 
ulated by hard gamma rays. In the remaining cases, the 
energy in this channel is not deposited in the atmosphere, 
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FIG. 1. Density distributions at ground level of e as a 
function of the distance to the shower axis. The primary 
zenith angle is 80° and E v = 10 20 eV. The error bars indicate 
the RMS fluctuations. 



and we will exclude such events from our consideration. 
The BH secondaries are injected with a primary zenith 
angle of 80° at 6.5 km above sea level (a.s.L), setting the 
observation level at 1.5 km a.s.l. All shower particles with 
energies above the following thresholds were tracked: 750 
keV for gammas, 900 keV for electrons and positrons, 10 
MeV for muons, 60 MeV for mesons and 120 MeV for 
nucleons. The geomagnetic field was set to reproduce 
that prevailing upon the Auger experiment. The results 
of these simulations were processed with the help of the 
aires analysis package. Secondary particles of different 
types in individual showers were sorted according to their 
distance R to the shower axis. We extract in a separate 
file all ^ and e^. 

In order to obtain a clear signature of BH production, 
one should be able to identify its subsequent cascade 
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in the whole cosmic ray sample. For large zenith an- 
gles (above 80°), an air shower initiated by a neutrino 
can be distinguished from that of an ordinary hadron 
by its shape. Ordinary hadrons interact high in the at- 
mosphere. As a consequence, at ground level the electro- 
magnetic part of the shower is totally extinguished (more 
than 6 equivalent vertical atmosphere were gone through) 
and only the muon channel survives. Besides, the shower 
front is extremely flat (radius > f 00 km) and the particle 
time spread is very narrow (At < 50 ns). Unlike hadrons, 
neutrinos may interact deeply in the atmosphere, trigger- 
ing showers in the volume of air immediately above the 
detector jllj]. The shower thus presents a curved front 
(radius of curvature of a few km), with particles well 
spread over time, O(ps). If primaries are mainly elec- 
tronic and muonic neutrinos (as expected from pion de- 
cays) two types of neutrino showers can be distinguished: 
"mixed" (with full energy) or "pure hadronic" (with re- 
duced energy), respectively. In the charged current inter- 
action of a v e , an ultra high energy electron is produced 
which initiates a large electromagnetic cascade parallel to 
the hadronic cascade. In contrast, the charged current in- 
teraction of a produces a muon which is not detectable 
at Auger. For the same total energy of the primary, the 
presence of a hard electromagnetic channel in BH pro- 
duction provides a clean signature when compared with 
the "pure hadronic" shower characterizing the inter- 
action. To analyze the differences between the BH-likc 
shower and ordinary v e shower, we mimic the latter as 
a superposition of a quark jet (equivalent to the set of 
hadrons listed in Table I) carrying around 20% of the 
original energy + a photon shower. Again, all particles 
in the sample are injected at 6.5 km a.s.l. and with a 
primary zenith angle of 80°. 

In Fig. 1 we show the e ± density at ground level (as 
a function of their distance to the shower axis), obtained 
from ordinary t^-shower and a BH-like shower. We set 
A/bh = 5 TeV, and E v = 10 20 eV. Then, for 4 < n < 7, 
from Eq.@ we get (N) w 5. At 50 m from the core, the 
ratio of the number of e ± in a BH-like shower to that in a 
typical v e shower is ~ 10~ 3 . At about 1 km from the core 
this ratio rises to ~ 10 -1 . Note that the differences far 
from the shower-core are also statistically significant for 
surface detector experiments like the Auger Observatory 
p6| . In Fig. 2 we show the resulting distributions of 
muons at ground level. This profile is seen to be a rather 
poor discriminator between BH and ordinary showers (of 
the same total energy), in spite of the fact that each 
of the four hadronic jets from black hole decay has the 
same energy as the single jet in the standard charged 
current interaction. There are sufficient muons produced 
by the lepton shower to largely close the gap between the 
profiles. 

In the presence of maximal f M /^ T -mixing, ^ T -showers 
must also be considered. However, since the mean flight 
distance ~ 50-E km/EeV, and the distance between posi- 
tion of first impact and ground is ~ 30 km, only r's with 
energy < 8 x 10 17 eV will decay. Thus, v T showers above 
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FIG. 2. Density distributions at ground level of ^ as a 
function of the distance to the shower axis. The primary 
zenith angle is 80° and E v = 10 20 eV. The error bars indicate 
the RMS fluctuations. 



this energy will be indistinguishable from showers. 

We comment briefly on the warped scenario. If the 
curvature spills into the extra dimensions the fundamen- 
tal Planck scale can be lowered all the way to ~ 1 TeV 
already for n = 1. For BH radii smaller than the cur- 
vature scale of the warp geometry, we expect validity of 
the flat space approximation. This leads to larger cross 
sections and jet multiplicities. The signature we have de- 
scribed remains robust when increasing the multiplicity 
up to about 15, because the properties of the electromag- 
netic shower depends almost exclusively on the ratio of 
hadronic to leptonic activity. 

In summary, cosmic neutrinos with horizontal inci- 
dence to the ground may interact with the earth at- 
mosphere producing BHs that decay instantaneously via 
Hawking evaporation. We have shown that such a reac- 
tion chain gives rise to deeply penetrating showers with 
an 'anomalous' electromagnetic component: about an or- 
der of magnitude bigger than ordinary j/^-showers and 
at least an order of magnitude smaller than ^.-showers. 
This represents a very clean signal. Our focus on a BH 
subsamplc in which the leptonic channel generates elec- 
tromagnetic showers will lower the event rate to about 
0.7/year (for n = 4, decreasing slowly for larger n). Thus, 
a 10-year collection of data at Auger could give signifi- 
cant statistics to test this phenomenon, yielding perhaps 
one of the early signatures of TeV-scale quantum gravity. 
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